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We study hadronization of the final state in a particle-antiparticle annihilation using a holographic 
gravity dual description of QCD. At the point of hadronization we match the events to a simple 
(Gaussian) energy distribution in the five dimensional theory. The final state multiplicities are then 
modelled by calculating the overlap between the Gaussian and a set of functions in the fifth dimension 
which represent each hadron. We compare our results to those measured in e + e~ collisions at LEP 
and PEP-PETRA. Hadron production numbers, which differ in range by four orders of magnitude, 
are reproduced to well within a factor of two. 



Introduction — Since the discovery of the AdS/CFT 
Correspondence [H, holographic gravitational theories 
have been studied to shed light on strongly coupled gauge 
theories. Phenomenological five dimensional (5D) models 
in this spirit (AdS/QCD) [2| also provide a quantitative 
description of the QCD meson spectrum which seem to 
work at the level of 10% accuracy. In this paper we will 
apply these tools towards another notoriously difficult 
QCD calculation: hadronization. 

Current Monte Carlo event generator models Q of 
hadronization are complex with many parameters which 
are tuned to data in some energy regime and have limited 
predictive power. A simpler understanding of the process 
would be a boon. Recently, progress has been made by 
assuming that after the quarks freeze into hadrons they 
may be described as a hadron gas in thermodynamical 
equilibrium [i| . Such models provide a surprisingly good 
description of the multiplicities of hadrons in jets across 
several orders of magnitude. 

Models of hadronization generically have two parts: 
predicting the initial yield of hadrons directly after anni- 
hilation and then allowing for decays of those particles in 
transit. Our model, like the thermal model [J], only ad- 
dresses the first part. Modelling the decays would involve 
a theory of branching ratios which we do not propose. In- 
stead we model this using the available branching ratio 
data from collider experiments. 

AdS/QCD models provide a weakly coupled 5D grav- 
itational theory that describes the hadrons made of con- 
fined quarks. Each hadron and its excited states (e.g. the 
"stack" of the p mesons: p(770), p(1450), p(1700) . . .) is 
described by a 5D field that shares the Lorentz properties 
and global symmetries of the hadron. In the gravitational 
theory one seeks solutions for these fields that separate 
the 3+1 dimensional dependence from the extra radial di- 
rection dependence, so they take the form g n {r)e lknX with 
k\ = M„. There are only regular solutions on the space 
for discrete masses corresponding to the meson masses. 
The functions g n (r) form an orthogonal basis. 

The dual necessarily only describes the regime of QCD 
where the gauge coupling is strong. When asymptotic 
freedom sets in at short distances and QCD is best de- 
scribed by free quarks and gluons the gravitational theory 



is expected to become strongly coupled itself and com- 
putation becomes impossible. Our description therefore 
is dependent on the distribution of energy at the match- 
ing energy scale where hadronization occurs. We will see 
that a very simple model reproduces the data very well. 

Hadronization Holographically — Our holographic 
model of hadronization will share with the thermal 
models the idea that at the point where the quarks 
form hadrons, the energy of the event is democratically 
available to all hadronic channels. We describe the 
initial condition as some deposition of energy into the 
5D model's stress-energy tensor. The radial dependence 
of the stress-energy dependence can be expanded in 
terms of the functions g n and will determine the relative 
multiplicities of each particle in a hadron stack. The x 
dependence will determine the energy and momentum 
of the hadrons. In this paper we will simply concentrate 
on the multiplicities. There should be some kinematic 
limit on the maximum mass of a state produced in the 
shower and in principle this could vary with the centre 
of mass energy of the event. In practice we shall simply 
include all known states with mass below 1.7 GeV (but 
excluding the ao(980) triplet - it is thought to be a 
bound kaon state). Above this value the experimental 
data on the full spectrum becomes patchy. In addition, 
the high mass states in this range are only produced 
with very small multiplicities and have a minimal effect 
on the lighter particle results. 

The result for the multiplicities depends on the choice 
of the function expanded in terms of the g n and this rep- 
resents the matching to the underlying asymptotically 
free QCD dynamics. The simple guess we will employ is 
that we should treat all hadronic channels equally and 
pick a Gaussian for the initial condition. The height of 
the Gaussian determines the absolute value of each par- 
ticle's multiplicity and hence is a free parameter which 
we fit (this parameter can be re-expressed as the aver- 
age energy per hadron, k, which we detail below). The 
width is also a free parameter which we fit. Finally the 
thermal model requires a suppression factor on the pro- 
duction of strange quarks. We too include a strangeness 
suppression factor 7 S which multiplies the Gaussian for 
each strange quark in the hadron (note that if mixing oc- 



FIG. 1: .The normalised holographic hadron basis functions for the 7r° stack (left) and the p meson stack (right) 



curs the strange quark content need not be an integer). 
7 S is another fit parameter in our analysis although we 
find it lies close to one. 

The final ingredient we require is a specific AdS/QCD 
model of the QCD hadrons. We will adapt a string theory 
derived model of chiral symmetry breaking that includes 
the vector mesons and pions That both the vector 
mesons and the pions are included is an important fea- 
ture: the mass spectra of the pseudo-Goldstone bosons 
are significantly different to all other hadrons in QCD. 
We will then assume that the g n functions associated 
with each hadron stack are not that different from the 
p-stack functions - we will simply reproduce them but 
with the mass of the lightest stack member tuned to the 
experimental value. Similarly the pion stack can be used 
to reproduce the towers of states associated with each 
pseudo-Goldstone of the chiral symmetry breaking (i.e. 
the pions, kaons and eta meson). The relative weighting 
is parameterised by R, a measure of gy M N, which we fit. 

Putting it all together, we compare our results to that 
detected at LEP (y/a = 91.2GeV) and PEP-PETRA 
(\/s = 29GeV), and find that we can reproduce yields 
that vary over four orders of magnitude to within a fac- 
tor of two. It is expected that the biggest source of error 
comes from the choice of holographic dual. 

Holographic Hadron Basis Functions — Our holo- 
graphic model of the hadron spectrum is based on a 
string theory construction of a QCD-like gauge theory 
consisting of a deformed D3 brane geometry with quarks 
included through probe D7 branes - the precise details 
do not need to be understood to follow this letter, but 
computational details are in the appendix and a more 
detailed analysis can be found in [5J. The model realizes 
dynamical chiral symmetry breaking. The mass of the 
lowest lying p meson can be dialed by choosing the con- 
formal symmetry breaking scale in the model, and the 
mass of the lowest lying pion can be dialed by choosing 
the asymptotic quark mass. The excited states in both 
stacks are then predicted: m p * = 1737 MeV, m„» = 1701 
MeV (c.f. experimental values of 1459 and 1300 MeV re- 
spectively). So whilst they don't precisely reproduce the 
experimental values the pattern is at least roughly right. 

As stated above, we assume that the g n functions asso- 
ciated with each hadron stack are not that different from 
the p-stack functions, and simply rescale the r coordinate 



such that the mass of the lowest member of each stack is 
correct. 

The functions representing the pions, henceforth de- 
noted as f n , are also rescaled. This time we dial the 
asymptotic quark mass such that the mass of the lowest 
member of each stack is correct. We plot the first three 
functions for f n , g n in figure [U 

The derivation of these functions is briefly reviewed 
in the appendix. We stress that here we seek to de- 
scribe hadronization not to provide a complete holo- 
graphic model of QCD. Hopefully the functions we use, 
based on the string model, are a reasonable phenomeno- 
logical basis. 

Overlap Computation of Multiplicities — With our 
holographic functions f n ,g n in place for each hadron 
stack we can now proceed with computing the expected 
initial yield in a hadronisation event. 

We assume that all the five dimensional fields ^(r) 
(e.g. each component of a gauge field A M describing the 
p mesons) have a common initial condition of a Gaussian 
centred at r — and of width A. To find the multiplicities 
of each stack member we compute 

/>oo 

Cn = / ^(r)w(r)g n (r)dr (1) 
Jo 

where w(r) is the weighting function associated with the 
basis functions g n . The Gaussian receives most support 
from the lowest mass states, whilst very highly excited 
states have no overlap with the Gaussian because of their 
highly oscillatory nature in the IR. 

The multiplicity is simply given by c„ multiplied by 
(2J+1) where J is the spin of the hadron. 

There are a number of special cases, which we address 
in turn. 

Pseudo-Goldstone Bosons — The pseudo-Goldstones 
are described by a separate holographic field, 9, in our 
model and are represented by functions f n (r). The rel- 
ative contributions the fields make to the stress-energy 
tensor are 

T rr ~ Ai(r)(M) a + A 2 (r)rr (cU M )(cU.) + . . . (2) 

Ai, A2 are functions of r, calculable from the holographic 
model. So that these fields see the same contribution to 



the stress energy tensor we therefore rescale the Gaus- 
sian. If our standard Gaussian is \&(r) then the pion 

field sees / dry/j*d r ¥. 

Strangeness Suppression Factor — Since the underlying 
asymptotically free dynamics may distinguish the strange 
quark from the up and down quarks, we also multiply the 
Gaussian by a factor of (j s ) (T where a is the strangeness 
content of the stack. j s is then the second fit parame- 
ter in our procedure. Note that this procedure is rather 
crude because different members of a stack may mix 
to varying degrees with other states. For example the 
t;(548) has 32% strangeness content while the n** has 
100%. In these cases we set 7 S by the strangeness con- 
tent of the lightest member of the stack. 

Height of Gaussian — The normalisation of the Gaus- 
sian tells us the relative multiplicities of the various 
hadrons in an event. An overall multiplicative factor n 
sets the absolute number of each species and we fit this 
value, k determines the total number of final state parti- 
cles (before allowing for decays in transit to the detector) , 
and hence we express it as the average hadron energy in 
the collision. 

A Fourth Parameter — Our choice of holographic dual 
also contains a free parameter, R, which sets the 't Hooft 
coupling in the gravity dual. We fit it to the data. How- 
ever, R is not in the same class as A,7 s ,k. R has a 
sound theoretical background, and would be predicted if 
the holographic dual to QCD was known. 

Decay in transit — Once we have calculated the ini- 
tial yield of hadrons, we then have to allow for decays of 
the particles in transit from the interactionpoint to the 
detector. Branching ratios are taken from [6], and parti- 
cles that can be detected at LEP (whose results we will 
compare to) are set as stable. All the other particles are 
allowed to decay through the decay channels until they 
reach one of the stable particles. In this way we get a list 
of numbers which is what our model predicts would be 
seen at LEP. 

Predictions — We compare our results both to e + e~ 
collisions performed at LEP (^/s = 91.2 GeV), and at 
PEP-PETRA (Vs = 29 GeV). Average multiplicities of 
various hadrons have been compiled in Q , and we repro- 
duce them here, along with our results, in table HI For 
both sets of results we have performed a four parameter 
fit so as to minimise the rms error. The pion predictions 
are expected to be low because our holographic model 
predicts over-massive excited states. Thus yields of par- 
ticles such as 7r(1300) and p(1450) are unnaturally sup- 
pressed which would otherwise be expected to give sig- 
nificant contributions to the pion multiplicities. Despite 
this, the fits are very good, with rms errors of 34% and 
46% for y/s = 91.2 GeV and ^ = 29 GeV respectively. 
The rj yields aren't ideal, but previous models have had 
the same problem Q. Furthermore, our model is ex- 
pected to have poor predictive power for the 77' meson: 



TABLE I: Results of the model for hadron yields at ^/s = 91.2 
GeV (centre column) and \fs = 29 GeV (right column). The 
relevant 4 free parameter values are shown in the final row. 
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0.03 


0.01 


0.02 
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0.172 


0.19 


0.13 


0.0983 
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0.0120 


0.0094 


0.0089 


0.0083 




0.012 


0.012 


0.0088 


0.0083 


= *0 


0.0040 


0.0033 


n/a 


n/a 


fl 


0.0011 


0.0014 


0.0008 


0.007 


A(MeV),K(GeV) 


150, 4.96 


152, 2.35 




0.97, 2.6 


0.97, 2.4 



technically it is a pseudo-Goldstone boson, but instan- 
ton effects cancel out this effect. Our model just treats it 
as non-Goldstone boson, which is probably over simplis- 
tic. In addition, it mixes heavily with 77(548): our model 
contains no good parameterization of mixing. 

The extra inaccuracy for the PEP-PETRA matching 
comes, in large part, from the f2 yield. This has a very 
large experimental error (see 0),' and so we may not be 
matching to the correct value: on comparing the PEP- 
PETRA and LEP experimental yields, a decrease by a 
factor of 10 adds doubt to the accuracy of the PEP- 
PETRA measurement. 

Conclusions — We assumed that every hadron in QCD 
can in principal be represented by a function in the r co- 
ordinate of the 5D holographic theory of QCD. We then 
proposed that hadronization can be modelled by hypoth- 
esising that the initial yield (that is before the particle 
created starts decaying) for any hadron is given by the 
square of the overlap between the function which repre- 
sents the hadron, and a Gaussian, centred at the origin, 
with a width of A. In addition we have two other pa- 
rameters in the theory; a strangeness suppression factor 
to account for the heaviness of the strange quark, and 
k which determines with what energy the particles leave 
the interaction point. 

With the full holographic dual to QCD currently un- 
known, we made some reasonable assumptions to achieve 
a full set of functions which represent every hadron. We 
then compared the results to e + e~ collisions made at 



LEP and PEP-PETRA. The results are surprisingly good 
suggesting the broad framework is correct. 

The model of hadronization presented here is applica- 
ble to all particle-antiparticle annihilation events, where 
the fireball after the collision has no residual quantum 
numbers. Broadening this model to include events such 
as deep inelastic proton-proton scattering, and heavy ion 
collisions would clearly be desirable. A natural way to 
do this would be to include enhancement factors on the 
multiplicities of stacks contributing to the quantum num- 
bers that are non zero in the final state. We leave such 
an analysis for the future. 

We would like to thank James Ettle and Francesco Be- 
cattini for their help in performing the branching ratio 
part of the calculation, and Ed Threlfall for his help in 
the preliminary stages of this paper. 

APPENDIX 

String Theory Progenitor — The phenomenological 
model used here is based on the AdS/CFT Correspon- 
dence realization of chiral symmetry breaking in 
That model consists of a dilaton flow deformed AdS ge- 
ometry 

b 4 ' 



ds 2 = H- 1 ' 2 flUx\ 



where H = f — 1, f 
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2<b 



6 4: 
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There are formally two free parameters, R and 6, since 
S = l/2b 4 and A 2 = 10 — S 2 . The parameter b sets the 
conformal symmetry breaking scale and we set it equal 
to one from this point onwards. We will use the p(770) 
mass to set an absolute mass scale. R sets the 't Hooft 
coupling which we treat as a free parameter. 

Quarks are introduced by including probe D7 branes 
into the geometry. We minimize the D7's world-volume 
in the spacetime around the D3 branes. This is encoded 
by the Dirac Born Infeld action in the Einstein frame of 
the D7 brane, which also contains a superpartner U(l) 
gauge field which describes the vector mesons. The full 
action is 



S 



which, expanded to second order gives 
S ~ [ d s £ e <#, V 3 5(l + ^ 2 )^ x 



1 



1 



-d a 6d a 6 



1 (2W) 



(1 + ct 2 ) 4(1 1 -• ' ® 

Our action is defined in 8D, rather than the usual 5D. 
But since we assume that the fields have no components 
in the 3-sphere (which is appropriate for duals to non- 
supersymmetric fields), J5]) can easily be recast as a 5D 
action. 

Substituting from the geometry above we can find the 
equation of motion for the radial separation, a, of the two 
branes in the 8 — 9 directions as a function of the radial 



coordinate r in the 4—7 directions. We first calculate the 
background solution as in [5] . The large r asymptotic so- 
lutions take the form o~q = m + c/r 2 with m representing 
the quark mass and c the quark condensate. The regular 
solutions have non-zero c even when m = and describe 
chiral symmetry breaking. 

Fluctuations of the brane about o~q in the angular di- 
rection, 8, of the 8 — 9 plane correspond to the pion. If we 
write 8(r,x) = 2ira / f(r) sin(fcx), the equation of motion 
for / is given by 

d ( e*G 



R 2 e*G 



<j 2 f +M 2 w f f = Q 



with wj — 



Hf 



1/4 



(r 2 



- 1 



(r 2 + ( x 2 ) 



2\2 



(G) 



(7) 



Finding regular solutions to (J6j) is simply a Sturm- 
Liouville eigenvalue problem, and hence we know that the 
regular solutions f n will form a basis under the weighting 
function wj. 

Similarly, if we write = g(r) sm(kx)e ll , the equation 
of motion for g is 



d 



g +M 2 w g g = (8) 



= R 2 e*gJ\ + a 2 Hf- 1 / i 



(9) 



Where once again g n form a basis under the weighting 
function w g . 

Stress-energy tensor — We also need to know the 
stress-energy contributions. 

T„ = — ?=JL (e*V=g£) (10) 



-gSg r 

We have to be careful, and use not ([5]), before 

expanding to second order. On doing so, we find 

A 2 g^g rr , , 

— = - — — (11) 



with Ai, A2 defined in equation ([2]). 



-det(P[ 3ah ]+27ra'e-^ /2 J F 1 Q6 ) * (4) {1] 
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